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We present an analytical solution for classical correlation, defined in terms of linear entropy, in 
an arbitrary d ® 2 system when the second subsystem is measured. We show that the optimal 
measurements used in the maximization of the classical correlation in terms of linear entropy, when 
used to calculate the quantum discord in terms of von Neumann entropy, result in a tight upper 
bound for arbitrary d ® 2 systems. This bound agrees with all known analytical results about 
quantum discord in terms of von Neumann entropy and, when comparing it with the numerical 
results for 10® two-qubit random density matrices, we obtain an average deviation of order 10“"*. 
Furthermore, our results give a way to calculate the quantum discord for arbitrary n-qubit GHZ 
and W states evolving under the action of the amplitude damping noisy channel. 

PACS numbers: 03.67.Mn,03.65.Ud 


Quantum entanglement plays important roles in many areas of quantum information processing, such as quantum 
teleportation and superdense coding [H-Q- Nevertheless, quantum entanglement is not the only form of quantum 
correlation that is useful for quantum information processing. Indeed, some separable states may also speed up 
certain quantum tasks, relative to their classical counterparts and many quantum tasks, such as quantum 

nonlocality iii and deterministic quantum computations with one qubit [lo| . can be carried out with forms of 
quantum correlation other than quantum entanglement. One such quantum correlation, called quantum discord, has 
received a great deal of attention recently (see 111 and references therein). Introduced by Ollivier and Zurek [1^ as 
the difference between the quantum mutual information and the maximal conditional mutual information obtained 
by local measurements 1^ [13|, quantum discord plays an important role in some quantum information processing 

Mill. 


Despite much effort by the scientific community, an analytical solution of quantum discord is still lacking even for 
two-qubit systems. Owing to the maximization involved in the calculation, there are only a few results on the analytical 
expression of quantum discord and only for very special states are exact solutions known. However, if instead of the 
von Neumann entropy one uses the linear entropy, the optimal measurements that maximize the conditional mutual 
information can be obtained analytically [2|. Here, we show that using these optimal measurements to determine the 
quantum discord in terms of the von Neumann entropy results in an excellent upper bound for the latter. Moreover, 
we show that this result gives a way to calculate the quantum discord for arbitrary n-qubit GHZ and W states, with 
each qubit subjected to the amplitude damping channel individually. 


RESULTS 

Classical correlation under linear entropy 

The usual quantum discord, in terms of von Neumann entropy, is defined as follows: let denote the density 
operator of a bipartite system composed of partitions A and B. Let and p^ = TrA{p^^) be 

the reduced density operators of A and B, respectively. The quantum mutual information, which is the information- 
theoretic measure of the total correlation, is defined as I{p^^) = S{p^) + S{p^) — S{p^^), where S{p) = —Tr(plog 2 p) 
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is the von Neumann entropy. Usually, the total correlation X is split into the quantum part Q and the classical part 
C, such that I = Q + C. The classical correlation of a bipartite state is defined as 

C(/^) = = nmx[5'(p^) - (1) 

i 

where the maximum is taken over all positive operator-valued measurements (POVM) {Pi\ performed on subsystem 
B, satisfying probability of i as an outcome, pi = Tr[(I^ (g) Pi)p^^{l^ ® P^)] where = 

Trs [(1^4 0 (I^ 0 )] 

is the conditional state of system A associated with outcome i, where is the identity 

operator on subsystem A. 

In this work, all POVM or projective measurements (PM) are taken on subsystem B. Finally, the quantum discord 
is defined as the difference between the total correlation and the classical correlation [13, [H: 

Qip^^) = I{p^^)-I^{p^^) 

= unnJ2p^S{p^'^)-S{p^\^), ( 2 ) 

i 

where S{p^^^) = S{p"^^) — S{p^) is the conditional entropy. 

To calculate our tight upper bound to quantum discord, instead of the von Neumann entropy one uses the linear 
entropy. The linear entropy of a state p is given by: 

S2{p) = 2[l-TTip% (3) 

In terms of the linear entropy ([3]), one can correspondingly define the conditional linear entropy, S 2 {A\B) = S 2 {p^^) — 
5*2 (p^), and the classical correlation 0 is written as: 

= max[S2ip^) - (4) 


where the measurements run over all POVMs Pi. 

Although the classical correlation and, consequently, the quantum discord ([2]) is extremely difficult to compute in 
terms of von Neumann entropy, the classical correlation (jT|) expressed in terms of linear entropy can be calculated 
analytically. In what follows we present the analytical formula for an arbitrary d 0 2 quantum systems. 

A qudit state can be written as p = (U-l-r • 7 )/d, where I^ denotes the dx d identity matrix, r is a (d^ — l)-dimensional 
real vector, 7 = ( 71 , 72 , ..., 7 ^ 2 - 1 )^ is the vector of generators of SU (d) and T stands for transpose. Consider a bipartite 
system, composed of a d-dimensional subsystem labeled A and a 2-dimensional subsystem labeled B. The bipartite 
state can be written as: 

P^^ = A 0 F(|yB,B)(VB,B|), (5) 

where \ Vb'b) is the symmetric two-qubit purification of the reduced density operator p® on an auxiliary qubit system 
B' and IF is the identity map on system B. Here, symmetric two-qubit purification means that the two reduced density 
matrices are equal, i.e. Vb' = Vb = p^ , and A is a a completely positive trace-preserving map which maps a qubit 
state B' to the qudit state A. Let cr = (cti, ( 72 , denote the vector of Pauli operators, r being a three-dimensional 
vector, |r| < 1 . As a qubit state can generally be written as p = (Id + r • o-)/d, the map A(p) is of the form 

A(p) = [Id + (Lr-Hs)- 7 ]/d, ( 6 ) 

where L is a {cP — 1) x 3 real matrix and s is a three-dimensional vector. L and s can be obtained from Eq. ([5]) and 
1 1 

Eq. dni). Let pB = be the spectral decomposition of pB- Then \Vb'b) = X] '/\\^i)\4’i) ^-nd A(|i)(jD, 

0 1=0 

f, j = 0,1, can be calculated by Eq. Therefore one gets A((Ti), i = 1,2,3, and the matrix Li^j = Tr(A((Tj).tTi). By 
the method used to calculate the classical correlation I^{p^^) of two-qubit states Q, we have: 

It{p^^)=\ma.{^^\^)S2{p^), (7) 

where Amoa:(L^L) stands for the largest eigenvalue of the matrix L^L. Eq. ([7]) gives the analytical formula for the 
classical correlation in terms of linear entropy for a general d 0 2 quantum state. Indeed, one only needs to find the 
eigenvalues of the matrix L^L. 
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Since, for a given state the reduced state \Vb'b) and the map A are fixed, the classical correlation can 
readily be computed in terms of linear entropy What concern us here are the optimal measurements that 

give rise to In fact, there is a one-to-one correspondence between all possible POVM measurements and 

all convex decompositions of p^ 0; namely, if p^ = ^ is the pure state decomposition of p'®, then the 

j=o 


following are the corresponding POVMs: 


Mo = {p^) =Po|'0o)(V'o|(p^) K (8) 

= (p^)~hil'ipi}('i/’il(p^)~K (9) 


where p^ is full-ranked. Otherwise, we can find the inverse of p^ in its range prthection and, from the optimal pure 
state decompositions of p^, we can get the corresponding optimal POVMs. In y, the authors have shown how to 
find the optimal decomposition of p^. First write p^ in its Bloch form: p^ = (/ -f tb • ct)/ 2. Let re -b Xj be the Bloch 
vector for the pure state decomposition \ijjj) of p^, where Xj = (VJ, Vj, VJ) and '^PjXj = 0, + Xj) = tb. 

j 3 

Hence, ||rB + ^j|| = l.Then = 1 — IItbIP — 2vB,^Xj — Without loss of generality, assume that 

is diagonal with diagonal elements > Xy > Eq. © becomes A 2 ,(l — HtbIP) + maxpj[(Ay — Xj;){X^)^ + 

(Az — Aa;)(VJ)^], which gets the maximum value when Xj = Xj = 0. There are exactly two solutions of the equation 

1 

||rB + Vjll = I. Hence the optimal decomposition of p^ reads: p^ = From the two pure states in the 

3=0 


optimal decomposition, we obtain the two optimal POVM measurement operators Mg and Mi. 

It is well known that to maximize the classical correlation it is necessary to use the most general POVM quantum 
measurement. As it is much more complicated to find the maximum in ([I} over all POVMs than over von Neumann 
measurements, almost all known analytical results are based on the latter. Indeed, only very few results are based 
on POVM [23,[2H. Here, we show that for the case of a bipartite qudit-qubit state, the classical correlation based on 
linear entropy is maximized over projective measurements (see proof in the appendix). This leads to our first theorem: 

Theorem 1. The classical correlation of a qudit-qubit state pab defined by running over all (arbitrary) POVM 
measurements is the same as the classical correlation defined by running over all projective measurements, i.e.. 


I. 


2 POVM 




Quantum discord under von Neumann entropy 


Theorem 1 implies that the optimal POVM in the classical correlation defined by Eq. ([4]) is in fact a projective 
measurement. This is very different from the case of classical correlation I'^{p^^) defined by von Neumann entropy, 
in which the classical correlation based on POVM could be larger than the one based on projective measurement 
2^ 211. This shows that, although von Neumann entropy and linear entropy have many properties in common, they 
behave quite differently in optimizing classical information. However, by using the optimal projective measurement for 
the classical correlation I^{p^^) based on linear entropy, we can get a tight lower bound for the classical correlation 
based on von Neumann entropy, and hence a tight upper bound for the quantum discord based on von Neumann 
entropy. This leads us to our second theorem: 

Theorem 2. The quantum discord based on von Neumann entropy has an upper bound: 


Q(/^) < I(/^) - [5(/) - ^p.5(/l*)]. 


( 10 ) 


where Pi = TT[{I^0Pi)p^^{I"^0P^)] is the probability of the measurement outcome i, = TrB[(/"^®Pi)p'^^(/"^® 
Pj)]/Pi is the conditional state of system A when the measurement outcome is i, and Pg s-nd Pi are the optimal 
projective measurement operators for of a given d®2 state pab- 

In fact, there is a connection between discord and entanglement of formation (EOF): the classical correlation 
I'^{p^^) can be obtained from EOF by the Koashi Winter Relation [23| . 

/-(/^) + P(p^c) = 5(/), (11) 

where is the original classical correlation of n® 2 state p"^^, E{pAc) is the EOF of state p"^^ , and p^'" is the 

purification of p"^^. It is important to note that, from theorem 2, we can get an upper bound of EOF for arbitrary 
rank two n® m state p^^. 
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Although the upper bound (IlOp of Q{p^^) is given by the optimal measurement of we show, by means 

of examples, that it is a surprisingly good estimate of Q{p"^^). 

Example 1. In Luo presented the analytic formula for the quantum discord Q{p^^) of the two-qubit Bell- 
diagonal state: p = (/ ® J 4 - ® o'i)/4. Let c = max{\ci\, |c 2 |, |c 3 |}. For this Bell-diagonal state, L = 

Diagjcf, c|, c|} and tb = 0. The two solutions of Xj are (1,0,0)’^ and (—1,0,0)^ when cf > max{c|,c|}: |V'o) = 
^{I+ai) and \ijji) = ( 0 , 1 , 0 )^ and ( 0 ,-1, 0 )^ when C 2 > max{cBC§}: IV'o) = ^il+cr 2 ) and \tjji) = ^{I-cr 2 ); 

(0,0,1)^ and (0, 0,—1)^ when c| >max{c 2 ,c|}: |^o) = and |'i/'i) = ^{I—as). It can be verified immediately 

that the optimal measurements for I^{p^^) are given by Mi = [I + CTfe)/2 and M 2 = (/ — crfe)/2, for c = Ck with 
k = {I, 2, 3}. It can easily be checked that our upper bound (ITOl) is exactly the same as the analytical results in ( 2 ^ . 

Example 2. In (j^ ? , |2^ the X-type two-qubit states are investigated: px = {I ® I + 2 : 3 (fra ® I) -\- y 3 {I ® (t^) + 
Si=i ® '2’i))/4j where 2 : 3 , j/ 3 , ti, t 2 and ta are defined such that p is a quantum state. It can easily be seen that 
our upper bound (fT0|) agrees perfectly with the analytical results obtained in [2^ (see Fig. I). 



FIG. 1: 


Quantum discord D = Q{px) for 2:3 = 0.1, 2/3 = 0.2, ti = 0.2, t 2 = 0.3. Here the results in 
numerical results and our upper bound in Eq. OT agree with high precision. 


m 


our 


Now, let us consider the following general two-qubit states, P 2 x 2 = 
and compare our analytical upper bound with numerical results. Fig.2 gives the quantum discord Q(p 2 x 2 ), for 
xi = 0.05, 2:2 = 0.1, 2:3 = 0.1, j/i = 0.15, 2/2 = 0.25, 2/3 = 0.2, ti = 0.2 and t 2 = 0.2 plotted against ^ 3 , such that p 2 x 2 
is a quantum state. Fig.3 shows the quantum discord Q(p 2 x 2 ) for X 2 = 0.1, 2:3 = 0.1, 2/1 = 0.15, 2/2 = 0.25, 2/3 = 0.2, 
ti = 0.2, t 2 = 0.2 and = —0.5, plotted against 2 : 1 , such that P 2 x 2 is a quantum state. It can be seen that our upper 
bound coincides very well with the numerical results. 



FIG. 2: Figure (a) shows quantum discord D = Q{p 2 x 2 )- Solid blue line shows numerical results and the red dotted 
line our upper bound. Figure (b) shows the difference between the numerical results and our upper bound. 


We have seen that the upper bound of quantum discord based on von Neumann entropy, obtained from the optimal 
measurements for the classical correlation based on linear entropy, is often exact. To test the precision of our upper 
bound generally, we calculated the difference between our analytical result and the numerical solution of quantum 
discord, for a set of 10® random density matrices of 2 0 2. In Fig.4, we plot the deviation A = QAnaiyticai — QNumericai 
against the number of occurrences. It can be seen that more than half of the randomly generated density matrices 
results in a precision greater than 10 “^, which demonstrates that our analytical result is a tight upper bound. 
Furthermore, in Fig.4, we show that more than 80% of the density matrices randomly generated lead to a precision 
greater than 10“®. Indeed, the percentage of density matrices with a deviation A greater than 6 x 10“® is less than 
0.1%. Here, in the horizontal coordinate of Fig. 4, 0..1 represents the interval from 0 to 1, and the same for 1..2, etc.. 
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FIG. 3: Figure (a) shows quantum discord D = Q{p 2 x 2 )- Solid blue line shows numerical results and the red dotted 
line our upper bound. Figure (b) shows the difference between the numerical results and our upper bound. 



FIG. 4: A as a function of number of occurrences for a set of 10® random 2 0 2 density matrices. 


Evolution of Quantum Discord under AD Channel 


Now we consider the evolution of quantum discord for arbitrary n-qubit GHZ and W states under an amplitude 

damping (AD) channel characterized by the Kraus operators Ei = ^ q q ^ i ?2 = ^ q ^ ■ We show 

that the related quantum discord based on von Neumann entropy can be analytically obtained from the upper bound 
given by Eq. m- 

First let us consider n-qubit GHZ states, with the first (n — 1) qubits subjected to AD channels individually. 
From Theorem 2, we get the optimal measurement operators {I 2 + o’z)/2 and {I 2 — o’ 2)/2 for classical correlation in 
terms of linear entropy, and the upper bound of quantum discord in terms of von Neumann entropy is then exact. 
Let Ml = U Pq and M 2 = U Pi U~^ be the two measurement operators, where Pq and Pi are the projective 
operators, U = tl 2 + yiOx + ?/ 20 'y + with = cos^ 6 and yl + y 2 = sin^ 0. Fig.5 shows that when 0 = 0 or 

9 = TT, has the minimal value, which coincides with the optimal measurement operators {I 2 + crz)l‘2 and 

{I 2 — crz )/2 for classical correlation based on linear entropy. 

For n-qubit W states with the first (n — 1) qubits subjected to individual AD channels, from Theorem 2 we have 
the optimal measurement operators (I 2 -b and {I 2 — o'a;)/2 or [I 2 + o’y)/2 and {I 2 — ijy)l2. The upper bound of 
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quantum discord obtained in terms of these measurement operators coincide with its lower bound in (^ . It follows 
that again we have the exact value of quantum discord ([2]). 

Alternatively, if the last qubit of an n-qubit W state is subjected to an AD channel, we have the optimal measurement 
operators {I 2 + (Jx)/‘^ and {I 2 — (Jx)/‘^ or {I 2 + o'y)/2 and {I 2 — o'y)/2, which also give rise to the exact value of discord 

m- 


CONCLUSIONS 

We have studied the quantum discord of qudit-qubit states. The analytical formula for classical correlation based 
on linear entropy has been explicitly derived, from which an analytical tight upper bound of quantum discord based on 
von Neumann entropy is obtained for arbitrary qudit-qubit states. The upper bound is found to be surprisingly good 
in the sense that it agrees very well with all known analytical results about quantum discord in terms of von Neumann 
entropy. Furthermore, for a set of 10® random density matrices, the maximum deviation found from the numerical 
solution was approximately 0.05 and the number of density matrices whose deviation was greater than 6 x 10“® was 
less than 0.1% of the whole set. Our analytical results could be used to investigate the roles played by quantum 
discord in quantum information processing. For classical correlation in terms of linear entropy, it has also been shown 
that the result for a qudit-qubit state, defined by running over all two-operator POVM measurements, is equivalent to 
that defined by running over all projective measurements. Furthermore, our results can be applied to investigate the 
evolution of quantum discord for arbitrary n-qubit GHZ and W states. Indeed, employing an important paradigmatic 
noisy channel, we present the quantum discord dynamics for the GHZ and W states when each qubit is subjected to 
independent amplitude damping channels. 
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Appendix 

[Proof of Theorem 1] Theorem 1 can be proved by using an approach similar to that used in [l[. It was proved 
by [3, for the classical correlation of ^ qudit-qubit state pab defined by running over all POVM 

measurements (here all POVM means we run over arbitrary POVM, that is, any measurement operators POVM), its 
optimal POVM measurement must be two operators POVM, that is, we can in fact restrict to two operators POVM. 
Then let M{x) and M{—x) be two optimal POVM operators, such that M{x) + M{—x) = I. Let Ho and Hi be the 
projective measurement operators. Ho -I- Hi = I. We can write the two POVM operators as M{x) = and 

M{—x) = where ao(a:) and ai{x) are the eigenvalues of M(x) and ao(—x) and ai(—x) are eigenvalues 

of M(—x). 

Given any qudit-qubit state pab, the POVM {M(x), M(—x)} performed on subsystem B will yield the post¬ 
measurement state PAjM^O- We have p(x)paimi^(x) = TtbIpabU 0 PA|nf ■ From the 


concavity property of linear entropy, we have the lower bound of the conditional linear entropy, 


P{y)S2iPA\M^{y)) 


y—x,—x 


= i: 

y—x,—x i—0 


y—x,—x i—0 


p{y) 


1 

y^P»»5'2(PAinf )■ 

i=0 


Thus, the conditional linear entropy derived from POVMs is greater than or equal to the conditional linear entropy 
derived from the projective measurements, on the all possible measurements basis. 

Let {nf} be the measurement basis that maximizes the classical correlation POVM operators 

{M[x),M{—x)}. Then, we have 


— ^2{p^) — '^2, P(2/)'^2(PA|MS(y)) 

y—x, — x 

< S2{p^) — y]]pi«S'2(p^|nf ) E ^2Pm(p"^^)> 

i 


since {nf} could be a non-optimal projective measurement of the classical correlation Hence, 

— -^ 2 Pm(p^^)! i-®-’ classical correlation under two POVM measurements is always smaller than 
or equal to the classical correlation under projective measurements. 

On the other hand, a projective measurement is a POVM. Hence, by definition, the classical correlation under 
two POVM measurements is always greater than or equal to the classical correlation under projective measurement. 
Therefore, we have proved that the classical correlation under two POVM measurements is equal to the classical 
correlation under projective measurement, I^poyyiiP^^) ~ ^ 
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